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ABSTRACT 

Using extensive HST imaging, combined with Chandra X-ray and Parkes radio data, we have detected the 
optical binary companion to a second millisecond pulsar (MSP) in the globular cluster 47 Tucanae. This faint 
(y = 22.3) blue (V -I = 0.7) star shows a large amplitude (60-70%) sinusoidal variation in both V and /. The 
period (3.19066 hr) and phase of the variation match those of the MSP 47 Tuc W (which does not have an 
accurate radio timing position) to within 0.5 seconds and L2 minutes respectively, well within the 1 a errors. 
The phase dependence of the intensity and color implies that heating of a tidally locked companion causes the 
observed variations. The eclipsing nature of this MSP in the radio, combined with the relatively large companion 
mass (> 0.13 M©) and the companion's position in the color-magnitude diagram, suggest that the companion is 
a main sequence star, a rare circumstance for an MSP companion. This system is likely to have had a complex 
evolution and represents an interesting case study in MSP irradiation of a close companion. We present evidence 
for another optical variable with similar properties to the companion of 47 Tuc W. This variable may also be an 
MSP companion, although no radio counterpart has yet been detected. 

Subject headings: binaries: general — globular clusters: individual (47 Tucanae) — techniques: photometric — 
pulsars: individual (PSR J0024-7204W) — pulsars: general 



dra and HST, have allowed the recent detection of two optical 
counterparts to cluster MSPs, PSR J0024-7203U, hereafter 47 
Tuc U, in 47 Tuc (Edmonds et al. 2001a, hereafter EGHOl), 
and PSR J1740-5340, hereafter NGC 6397 A, in NGC 6397 
(Ferraro et al. 2001). The MSP 47 Tuc U has a -0.15 Mq 
He white dwarf (WD) companion, a 10.3 hr binary period and 
small amplitude (0.004 mag in V) orbital variations caused by 
heating of one side of the companion (EGHOl). The '^0.45 
M0 companion in NGC 6397 A is either a subgiant or a heated 
main sequence (MS) star (Ferraro et al. 2001). The 32.5 hr, 
0. 1 mag orbital variations are caused by tidal distortion of the 
secondary by the NS, and are roughly sinusoidal with a (32.5/2) 
hr period. Another notable feature of NGC 6397 A is that the 
spectrum of the X-ray counterpart is relatively hard, and is sug- 
gestive of non-thermal emission. By contrast, 8 out of 9 of the 
47 Tuc MSPs which are bright enough to have useful spectral 
information are soft (Grindlay et al. 2002). 

Currently, NGC 6397 A appears to be unique among known 
MSPs in having a likely non-degenerate companion. The MSPs 
in 47 Tuc, for example, are either single, have likely He WD 
companions (e.g. 47 Tuc U), or have secondaries with masses 
of -0.02-0.03 M0 (CLFOO) which cannot be normal MS stars. 
As argued by Burderi et al. (2002) and Ferraro et al. (2001), the 
subgiant companion to NGC 6397 A may have been the star that 
recently spun up the MSP, consistent with the relatively small 
(characteristic) age for the system (D'Amico et al. 2001). The 
position of the star in the color-magnitude diagram (CMD) is 
roughly consistent with that of a mass-losing subgiant like that 
in the CV AKO 9 in 47 Tuc (Albrow et al. 2001). The equally 
interesting possibility presented by Ferraro et al. (2001) is that 
an MS star has been captured in an exchange interaction with 
an NS already spun-up to millisecond periods. Further work is 

' Based on observations with the NASA/ESA Hubble Space Telescope obtained at STScI, which is operated by AURA, Inc. under NASA contract NAS 5-26555. 
^ Harvard-Smithsonian Center for Astrophysics, 60 Garden St, Cambridge, MA 02138; pedmonds@cfa.harvard.edu; cheinlce@cfa.harvard.edu; josh@cfa.harvard.edu 
^ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218; gillil@stsci.edu 

'* Columbia Astrophysics Laboratory, Columbia University, 550 West 120th Street, New York, NY 10027; fernando@astro.columbia.edu 
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1. INTRODUCTION 

More than 50 millisecond pulsars (MSPs), half of all known, 
have been found in radio surveys of globular clusters, where 
the high stellar densities and interaction rates cause large num- 
bers of neutron stars (NSs) to be spun up to millisecond peri- 
ods (via binary production and subsequent accretion onto the 
compact object). Studies of these objects therefore offer in- 
sight into the formation and evolution of NS binaries (see, e.g., 
Rasio, Pfahl, & Rappaport 2000) and the frequency of stellar 
interactions in the dense cores of clusters, sampling different 
evolutionary channels compared to the disk population. 

Recent improvements in the sensitivity of the Parkes radio 
telescope have resulted in a dramatic increase in the number 
of cluster MSPs detected. In particular, the 1 1 known MSPs 
in the cluster 47 Tuc (Manchester et al. 1991; Robinson et al. 
1995) have recently been increased to 20 (Camilo et al. 2000, 
hereafter CLFOO), with 16 of these currently having accurate 
(milli-arcsecond) timing positions (Freire et al. 2001a; Freire 
2001). These enable studies of the cluster's gravitational poten- 
tial well (Freire et al. 2001a) and of the intra-cluster medium 
(Freire et al. 2001b). 

Another significant observational advance has been the abil- 
ity of Chandra to detect cluster MSPs in large numbers (Grind- 
lay et al. 2001a, 2002; D'Amico et al. 2002). Accurate (0."1- 
0"2) X-ray positions for the MSPs combined with the detection 
of large numbers of cataclysmic variables (CVs) and active bi- 
naries with both Chandra and //5r (Grindlay et al. 2001a and 
Edmonds et al. 2002a, in preparation) has allowed the radio. 
X-ray and optical data to be placed on a common astrometric 
frame, good to <0"1. 

These advances, and the superb spatial resolution of Chan- 
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needed to test whether the position of the star in the CMD is 
consistent with estimates of heating of the companion by the 
MSP (Ferraro et al. 2001). The relatively large amplitude of 
the ellipsoidal variations in NGC 6397 A suggests a high incli- 
nation orbit (Ferraro et al. 2001; Orosz & van Kerkwijk 2002), 
consistent with the detection of echpses in the radio (D'Amico 
et al. 2001). 

Five out of eight eclipsing MSPs known have very low com- 
panion masses of ~0. 02-0.03 M0, and short orbital periods 
(these include the 47 Tuc MSPs J, O, and R). The three ex- 
ceptions are NGC 6397 A; PSR B1744-24A in the cluster 
Terzan 5, with an orbital period of ~1.8 hr, companion mass 
^0.10 Mq, and displaying irregular eclipses (Lyne et al. 1990); 
and 47 Tuc W, an eclipsing MSP in 47 Tuc with a ~3.2 hr pe- 
riod but with a companion mass of ~0.15 Mq (CLFOO; see 
§2.1). 

We present here extensive //ST evidence for an optical coun- 
terpart to 47 Tuc W based on the detection of a faint, large am- 
plitude, optical variable with an orbital period and phase which 
perfectly match those of 47 Tuc W within the uncertainties. The 
X-ray counterpart detected with Chandra has properties which 
are similar to those of NGC 6397 A. This detection represents 
the best evidence yet obtained for an MSP with an MS compan- 
ion, given uncertainties about the evolutionary status of NGC 
6397 A. Following a brief summary of the radio data available 
on 47 Tuc W, we present the astrometry, absolute photometry 
and time series photometry for the 47 Tuc W optical companion 
in § 2, where we also present evidence for a 2nd possible MSP 
companion with similar optical and X-ray properties. In § 3 we 
discuss the impUcations of these results. 

2. OBSERVATIONS, DATA ANALYSIS, AND RESULTS 

2.1. Radio Data 

The 2.35 ms pulsar 47 Tuc W was discovered on 1999 Febru- 
ary 5 (MJD 51214) during a ~4 hr observation with the Parkes 
radio telescope (CLFOO). It has never been detected again 
(which happens for a few very weak pulsars that are detected 
seldom and only due to the focusing effects of interstellar scin- 
tillation), and hence has no position determined from radio tim- 
ing measurements. However, that one observation was suffi- 
cient to determine that it is part of a binary system with a period 
of ~3.2 hr, owing to the varying apparent spin period caused by 
Doppler shift. The times-of-arrival (TOAs) determined from 
the discovery data can be used to obtain precise orbital param- 
eters. 

We have used 40 good TOAs together with the TEMPO^ tim- 
ing software to measure an accurate barycentric pulsar spin pe- 
riod, projected semi-major axis of 0.2434(1) light-sec, binary 
period Ph = 0. 13295(4) days, and time of pulsar ascending node 
7b = 51214.284829(3) MJD, where all uncertainties given are 
conservative (lOx TEMPO), accounting for covariances in fit- 
ted parameters. The orbit is circular, with eccentricity < 0.001. 
With orbital phase i(j)b) zero defined at 7b, the pulsar is clearly 
eclipsed during 0.1 < (ph <0.4 (hereafter we redefine the fiducial 
point of orbital phase such that (ph = at T^ax = Tq+ [3/4]Ph, 
retaining the binary phase convention used by EGHOl). With 
the measured orbital parameters, and assuming an NS mass of 
1.4 M0, the pulsar mass function implies a companion mass 
> 0.13 Mq, with a mass of ~0.15 Mq for an orbital inclination 
of 60°. 

^ See http://pulsar.princeton.edu/tempo. 



2.2. Summary of Optical Data 

The optical data used in this paper consists of four different 
programs, GO-8267, GO-7503, GO-8219, and SM3/ACS- 
9028, summarized in Table 1. The program GO-8267 (PI R. 
Gilliland) involved an extensive set of 160 s Wide Field Plan- 
etary Camera 2 (WFPC2) exposures in F555W (636 images) 
and F814W (653) obtained over 8.3 days in 1999 July. A lim- 
ited number (28) of F336W exposures were also obtained (not 
shown in Table 1). The use of these data to search for tran- 
siting planets and to study binaries are described in GilUland 
et al. (2000) and Albrow et al. (2001) respectively. The pro- 
gram GO-7503 (PI G. Meylan) involved a short set of dithered 
F555W exposures with WFPC2 taken over three orbits about 
110 days after GO-8267. The GO-8219 (PI C. Knigge) data 
consisted of six blocks of eight 30 s exposures each using the 
Space Telescope Imaging Spectrograph (STIS) in CCD mode 
with the so-called CLEAR filter (two exposures of one visit 
were lost resulting in 46 total exposures). The high through- 
put of STIS/CCD resulted in comparable count levels for faint 
stars in 30 s compared to the 160 s F555W and F814W WFPC2 
exposures. The first visit of GO-8219 was about 65 days after 
GO-8267, while the remaining five visits were over a 10 day 
period some 401 days after GO-8267. Finally, SM3/ACS-9028 
consisted of 20 cosmic ray split observations of 1 20 s with the 
F475W filter using the High Resolution Camera (HRC) of the 
newly installed Advanced Camera for Surveys (ACS) on HST. 
These data were acquired with 20 large-scale dithers over 4 or- 
bits of HST as caUbration of the camera's geometric distortion, 
and were obtained in 2002 April about 1002 days after GO- 
8267. 



2.3. W29 and Chandra-HST Astrometry 

The astrometry between Chandra and HST, based on the GO- 
8267 data, has already been described in EGHOl and Edmonds 
et al. (2002b). A partial set of optical identifications was de- 
scribed by Grindlay et al. (200 1 a) and the full set is described in 
Edmonds et al. (2002a, in preparation). This paper focuses on 
the optical identification for the X-ray source W29 from Grind- 
lay et al. (2001a). We first establish the similarity of this source 
to some known X-ray counterparts of MSPs, and then present 
the Chandra-HST asXxomeXxy. 

The low luminosity and relatively hard X-ray spectrum of 
W29 are similar to those of U12 (the X-ray counterpart of NGC 
6397 A; Grindlay et al. 2001b) and 47 Tuc J, both eclipsing 
MSP binaries (see Grindlay et al. 2002). The X-ray spectrum of 
W29 can be best fit with a power law (photon index = 1 .8±0.6 
when fixing Nh to the nominal cluster value; xj- = 1-55) and 
the luminosity in the 0.5-2.5 keV band is 7.8 x lO-'^erg s"'. 
By comparison, for U12 a power law fit gives a photon in- 
dex of 1.6±0.3 (when fixing Nh to the cluster value) and a 
luminosity of 4 x lO'^'erg s~'. Because both W29 and U12 
are relatively faint sources (41 and 66 counts respectively), 
these X-ray sources are also consistent with ^6 keV thermal 
bremsstrahlung spectra (x^ = 1-94 for W29). The X-ray coun- 
terpart of 47 Tuc J, W63, has only 10 counts (L^ = 1 .8 x 10'°erg 
s"') and is too faint for useful spectral fitting, but has hardness 
ratios consistent with a power law of photon index between 1 



An MSP Optical Counterpart in 47 Tuc 



3 



and 1.5 (Grindlay et al. 2002). 




Fig. 1.— HST finding charts for W29opt, using tlie GO-8267 (V) and 
SM3/ACS-9028 (F475W) data. Tlie mean direct image from GO-8267 
(WFPC2) is sliown in Fig. la and the mean difference image is shown in 
Fig. lb (see text for definitions). The ACS/HRC image for W29opt near max- 
imum ((j>i, = 0.0±0.16) is shown in Fig. Ic and the corresponding image near 
minimum (0j = 0.5±0. 16) in Fig. Id. The 5 a (Of'42) error circle for W29 is 
shown, as is the scale and orientation of the HST images. The variable, optical 
counterpart for W29 (denoted W29opt) is clearly visible near the center of the 
X-ray error circle in the mean difference and maximum ACS images (labeled 
with an arrow), but is barely visible in the mean direct and minimum ACS im- 
ages. The initial identification for the ^^3.2 hr variable (PC1-V46) is labeled 
with 'X'. 

A plausible optical counterpart for W29 was identified as a 
nearby 3.187 hr period variable. The variations in this very 
crowded region were initially believed to arise from a relatively 
bright star {U = 18.37, V = 17.79 and / = 17.13), identified as 
a non-eclipsing contact binary with an orbital period of 6.374 
hr, and an amplitude of ^^0.002 mag (PC1-V46 from Albrow et 
al. 2001). This star is located 0."25 (3.0 <t) from W29, where 
the astrometric uncertainty standard deviation is a combina- 
tion of positional errors for W29 and small systematic errors in 
transforming between the Chandra and //ST coordinate frames. 
However, a careful examination of the directly observed images 
('direct images'; a term used to describe both individual expo- 
sures and mean images) and difference images^ showed that 
the bright star identified as PC1-V46 was not the true variable. 
Instead, there is a faint, nearby star (only 0"19 away) with a 
much larger intrinsic variability amplitude. This star was not 
included in our original star lists because it is faint and very 
crowded (it lies only 3."8, or 16% of the core radius, r^-, away 
from the cluster center), but it fell within, and hence contam- 
inated, the aperture used for PC1-V46. The new variable will 
be designated W29opt since it is only 0."07 (0.8 a) away from 
W29. The likelihood that this positional match is a coincidence 
is very low. On the PC chip there are 39 variables from Albrow 
et al. (2001) with V fainter than the MS turnoff, so the proba- 

^ difference images are generated by subtracting from each individual direct ima] 
Albrow et al. 2001 



bility of a chance match between W29 and an optical variable 
is only ^ 5 x 10"^. The position of W29opt in the astrometric 
coordinate system of the MSPs (based on the JPL DE200 plan- 
etary ephemeris; Freire et al. 2001a) is given in Table 2. 

The time series for W29opt is approximately sinusoidal with 
an amplitude of 70% in the V-band and 60% in the /-band. 
More details on the time series will be given in § 2.5. Here, 
we use the phase information for W29opt to construct mean di- 
rect and difference images for this variable. Figure la shows 
the average of direct images within ±0.07 in orbital phase of 
cj)h = 0.0 (maximum light) and within ±0.07 of (pb = 0.5 (min- 
imum light) in the W29opt time series, where we use the bi- 
nary phase convention given in § 2.1. This provides an average 
intensity direct image for the variable without over-sampling. 
The variable W29opt is barely visible about 1 pixel below the 
center of the 5 a X-ray error circle for W29. Clearly, W29opt 
is located in a crowded region of 47 Tuc even at HST resolu- 
tion. Although there are no other identifiable stars within a 0''1 
radius of W29opt, within a radius of 0"2 there are three stars 
with a total intensity in F555W over 200 times that of W29opt, 
and within a 0."5 radius a total of 13 stars have a total intensity 
nearly 800 times that of W29opt. Since W29opt has a color com- 
parable to cluster turnoff stars (see § 2.4), while most stars are 
redder than this, the relative brightness of nearby stars is some- 
what higher yet in F814W or the STIS CLEAR filter bandpass. 

Figure lb shows the sum of all difference images within 
±0.07 of (pb = 0.0, minus the sum of all difference images 
within ±0.07 of (j>b = 0.5, giving the mean difference image. 
Non-variables disappear except for noise. Note the clear de- 
tection of W29opt near the center of the X-ray error circle and 
coincident with a faint, relatively blue star just visible in the 
mean direct image. The large variability of W29opt is shown by 
comparing the two ACS/HRC images shown near maximum 
(c/)b = 0.0±0.16) and minimum ((j)b = 0.5±0.16). The absolute 
photometry for this star will now be described. 

2.4. Photometry 

Because of the extreme crowding described in § 2.3, the ab- 
solute photometry for W29opt requires special attention and we 
have used two different methods to estimate the U, V and / 
magnitudes for this star: 

(/) We constructed direct images for W29opt at maximum {(pb 
= 0.0±0. 16) and minimum {(pb = 0.5±0. 16) in V and /, and sub- 
tracted the minimum image away from both the maximum and 
the mean images. In each of these two residual images, W29opt 
was much more clearly visible than in any of the direct images, 
and crowding was no longer a limiting factor in the photom- 
etry. After carrying out point spread function (PSF) fitting to 
the W29opt neighbors in the mean image, we applied this PSF- 
fitting to W29opt in the residual images, and used the measured 
light curves to derive V and / magnitudes for W29opt at (pb = 0.0, 
(pb = 0.5, and mean values averaged over the complete GO-8267 
dataset. For U a mean value only was derived using this tech- 
nique, because of the much lower signal-to-noise ratio. Mean 
values of f/ = 23.8 ± 0.5, V = 22.3 ± 0.2 and / = 21 .6 ± 0.3 were 
derived (see Table 2), where the errors are a combination of un- 
certainties from the PSF-fitting and estimated absolute errors in 
the conversion from HST to Johnson magnitudes. 

(ii) The second technique involved use of the deep, oversam- 
pled images, the accurately known position for W29opt, and 

;e the mean over-sampled image evaluated at the appropriate dither position; see 
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local sky determinations, to estimate the peak intensity levels 
for W29opt in the three different filters. We then scaled these 
values to Johnson UVI magnitudes using nearby stars and the 
PSF-fitting results of Gilliland et al. (2000). Mean values of 
f/ = 23.7 ± 0.7, y = 22.4 ± 0.2 and / = 21 .7 ± 0.3 were derived, 
in good agreement with the values from (;). 




-2 2 -0.5 0.5 1 1.5 

U-V V-I 



Fig. 2.— HST CMDs for GO-8267 ba.sed on PSF-fitting. The variable 
W29opt and the counterpart to the MSP 47 Tuc U ('U'; EGHOl) are labeled. 
The thicker lines show Bergeron, Wesemael, & Beauchamp (1995) models 
for 0.5 Mq (dashed line) and 0.6 Mq (solid line) CO WDs. The thinner 
lines are SereneUi et al. (2001) models, for He WD masses of 0.406, 0.360, 
0.327, 0.292, 0.242, 0.196 Mq, with mass decreasing towards the red, and 
0.169 M0 (dot-dashed line). Cooling ages are as shown. These models have 
been plotted assuming that (m-M)o = 13.27, the mean of the 9 distance mod- 
ulus values reported in Zoccali et al. (2001), E(B-V) = 0.055 (Zoccali et 
al. 2001), Au/Av = 1.51 and Ai/Ay = 0.60 (Holtzman et al. 1995). Using 
(m-M)o = 13.27 the distance to 47 Tuc is 4.5 kpc. 

Figure 2 shows the V vs fZ-V and V vs V-/ CMDs for 
the PC chip, with W29opt labeled. The known counterpart to 
MSP 47 Tuc U (EGHOl) is also shown. All of the points in the 
CMDs, including W29opt (where the photometry from [;] was 
used), had magnitudes which were self-consistently determined 
by PSF fitting. The various plotted Unes are the CO WD cooling 
tracks of Bergeron et al. (1995) and the He WD cooling tracks 
of SereneUi et al. (2001). Note that the colors of W29opt ap- 
pear to be inconsistent with those of any He WD cooUng track, 
unUke the case for 47 Tuc U. 

2.5. Time Series 

2.5.1. GO-8267 data 

The original time series analysis, upon which the initial pos- 
sible match of W29opt with the MSP 47 Tuc W was made, used 
only the GO-8267 data. Here, we applied the standard differ- 
ence image technique (see Gilliland et al. 2000 and Albrow 
et al. 2001) using a circular aperture, for this crowded target 
star, with an area of only 5 pixels (1 pixel = 0"045 x 0."045) 
to minimize contamination from neighboring companions. The 
top panel of Figure 3 shows the full time series for W29opt, 
and the middle and bottom panels show the GO-8267 data with 
a sinusoidal fit superimposed (the /-band Ught curve has been 
scaled by a factor of 1.17 so that V and / have the same am- 
phtude). The times displayed in the x-axis were derived by 
extracting the Modified JuUan Date (MJD) of the observations. 



given in Universal Time, from the HST¥YYS headers, and con- 
verting to Heliocentric Julian Date (HJD) using the IRAF task 
rvcorrect^. The y-axis units are fractional intensity, defined 
as (intensity— <intensity>)/<intensity> (where <intensity> is 
the mean intensity). This system is used because it is linear (un- 
hke the magnitude system) and gives light curves for W29opt 
which correctly appear sinusoidal for an intrinsic sinusoidal 
modulation. 
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Fig. 3. — ifSTtime series for W29opt. The top panel shows the fuU set of 
time series analyzed here (the different datasets are labeled), extending from 
1999 July to 2002 April. The times shown are extracted from the HSrFITS 
headers. The middle panel shows the GO-8267 data and the bottom panel 
shows a close-up of a typical sample of the GO-8267 data (a sinusoidal fit to 
the data is shown in each case). The /-band data has been scaled so that its am- 
plitude equals that of the V-band data. In the bottom panel different symbols 
have been used for the two filters. 

We then performed a simultaneous least-squares fit of the 
GO-8267 y and / data to a sinusoid, allowing the amplitude, pe- 
riod and time of optical maximum ijrma = T[<j)b = 0.0]) to vary. 
The following values were obtained: amplitude = 0.714(9), 
n = 0.132919(16) d, and 7;„^ = 51366.87064(28) MJD (all 7;„^ 
values are given in barycentric MJD to within a small error of 
1-2 s). Because of the difficult photometry for W29opt there 
is an additional ~20% error for the amplitude resulting from 
assumed normalization using the mean absolute magnitude. 

We now compare these results with the measured radio pa- 
rameters of 47 Tuc W. The orbital period of 47 Tuc W (see §2.1 
and Table 2) differs from the period of W29opt by only 0.00003 
d (2.7 s, a 0.02% relative difference), a difference of 0.8 a when 
using the larger (radio) error. The radio value for T„ax (when the 
heated hemisphere of W29opt is facing towards the Sun and is 
expected to be brightest, in a model where W29opt is heated by 
the pulsar wind) is 51214.38454(3) barycentric MJD (Table 2). 
By projecting the optically-derived T,„a.K backwards in time by 
1147 orbits, we derive W = 51214.412(18) MJD. This differs 
from the radio-inferred T^ax by (0.028 ± 0.018) d, a 40 minute 
difference (0.2 in phase). These results suggest that W29opt 
may be the companion to 47 Tuc W, although the phase errors 
are still relatively large. These errors are significantly reduced 
in § 2.5.2 with the use of archival data. 

The phase plots and Lomb-Scargle (Scargle 1982) power 
spectra for W29opt in the V and / bands are shown in Figure 4. 



' HJD differs from barycentric MJD by a variable heliocenter-barycenter light travel time of ~l-2 s, in addition to a trivial 0.5 d offset due to a different zero-point 
definition of the two time scales. 
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Fig. 4. — Phase plots and power spectra for W29opt, using the GO-8267 data. The left side shows the V-band data and the right side the /-band data. The period 
used to generate the phase plots is slightly (0.02%) higher than the value derived from the GO-8267 data alone, because of improvements based on archival analysis 
(described in § 2.5.2). The squares show the mean fractional intensity in 25 phase bins between 4>i, = 0.0 and 4>i, = 1.0, where the error bars give the 3 ct errors (cr is 
the standard error in the mean). The solid line shows a sinusoid with the fitted amphtude, period and phase. Note the differences between the V-band data and the 
sinusoidal model near (pi, = 0.0 and (ph = 0.5 (see § 2.5.3). The vertical dotted line in the power spectrum shows the HST ovhital period (96.4 min) and the horizontal 
dashed line shows the Lomb-Scargle power corresponding to a false-alarm probability (Scargle 1982) of 1 x 10^. The power spectrum peaks near the main period 
are aliases caused by the /fSr window function. 



Note the similarity of the time series (for both filters) to a sinu- 
soid. The power spectrum shows a large peak at the frequency 
value corresponding to the fitted period, with no evidence for a 
peak at twice the fitted period. The modulation is thus consis- 
tent with arising from a single heated face of the secondary in a 
circular orbit. 

2.5.2. Use of archival data 

With the above solution as a starting point, we used images 
from the archive (GO-7503, GO-8219 and SM3/ACS-9028) to 
refine the orbital solution for the period and phase of W29opt. 
The errors for the period and phase of W29opt from GO-8267 
are small enough that the archival data can be included in the or- 
bital solution without any cycle count ambiguity. For example, 
the 401 day separation between GO-8267 and GO-8219 gives a 
phase uncertainty of 0.048 d (0.36 of one orbit), when using the 
error quoted above for the period. Therefore, there is potential 
to significantly improve our orbital solution using these extra 
data. 

The small set of F555W data from GO-7503 was taken in one 
contiguous period and spans about 0.85 in phase for W29opt. 
The combination of multi-pixel dithers for the small number of 
exposures from GO-7503, coupled with the strong variability of 
W29opt, implied that the GO-7503 data itself was not sufficient 
to generate difference images relative to an internally generated 
over-sampled mean image. To obtain difference images neces- 
sary for extracting variations of W29opt from the GO-7503 data 
we used the mean, over-sampled image from GO-8267. Posi- 
tions were fitted to about 50 stars within 7" of W29opt in both 
the over-sampled reference image and each individual exposure 
from GO-7503. Transformation coefficients were derived in- 
cluding x,y offsets, rotation and independent plate scales in x,y 



from the two sets of star positions. Then cubic-spline interpola- 
tion in the (essentially noiseless) over-sampled reference image 
from GO-8267 was used to evaluate an image at the position 
corresponding to each GO-7503 exposure. Difference images 
were then formed and the photometry at W29opt was obtained 
as a 5 pixel sum in each difference image at the calculated po- 
sition of W29opt. Since this is a very crowded region, and the 
two data sets had an orientation difference of 113°, there was 
a net (arbitrary) zero point offset removed to force consistency 
between the relative photometry of the GO-7503 and GO-8267 
datasets. The same scale factor (mean intensity of W29opt in the 
direct image) was used for GO-7503 (compensating for relative 
exposure times) as for GO-8267, taken with the same filter and 
detector. 

As with the GO-7503 data, there were not enough indepen- 
dent dithers and exposures for GO-8219 to generate good dif- 
ference images from these data alone. The GO-8219 data were 
rotated by ~54° relative to GO-8267, and were taken with a 
different instrument/filter combination. Again, fits were per- 
formed for about 50 stars within 7" of W29opt in both the PCI 
over-sampled reference image and the individual STIS expo- 
sures. Since the STIS CLEAR filter has high response spanning 
the WFPC2 F555W and F814W filters, we evaluated reference 
images using cubic spline interpolation from each of these two 
over-sampled images. The match to STIS data was then eval- 
uated as a linear combination of F555W and F814W reference 
images and zero point offset, providing a best fit to 50 fiducial 
stars and several nominal regions representing sky. A mean 
contribution of 0.437 times V and 0.590 times the / images 
from PCI were adopted to best represent the STIS CLEAR data 
and difference images were formed. To further match the STIS 
and WFPC2 data images we evaluated and appUed convolution 
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kernels attempting to match the relative PSFs. A best solution 
was found using a convolution kernel that broadened individual 
STIS PSFs by 10-15% relative to those for the PCI references. 
The new difference images were then used in a standard way to 
extract relative intensity values for W29opt using 5 pixel sums 
at the calculated position of W29opt in each frame. Here, rela- 
tive to the WFPC2 data sets, there is an additional uncertainty 
in normalization of the relative photometry for W29opt at the 
~25% level. 



(a) GO-7503: original solution 
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Fig. 5.— Phase plots for archival data from GO-7503, GO-8219 and 
SM3/ACS-9028. Figure 5a shows the sinusoid derived from the GO-8267 data 
only (the 'original solution') applied to GO-7503, and Figure 5b shows the re- 
vised solution based on the use of all four datasets. Similarly, Figures 5c and 
5d show the original and revised solution for GO-8219 (where 'CLEAR (1)' 
refers to the 1999 data and 'CLEAR(2)' refers to the 2000 data) and Figures 
5e and 5f the results for SM3/ACS-9028. Note the dramatic improvement in 
the results for GO-8219 and SM3/ACS-9028 based on the revised solution. 

For SM3/ACS-9028 the improved spatial resolution and PSF 
for HRC/ACS compared to WFPC2 permitted direct aperture 
photometry for W29opt, rather than difference image analysis 
(the star is clearly visible in the direct images). The data were 
dithered by close to half the detector size as part of this calibra- 
tion program, but W29opt was found in all of the frames. 

We then applied the 'original' sinusoidal solution from the 
8.3 d of GO-8267 data to the archival time series. Plots of the 
GO-7503 data (-110 d after GO-8267) using the above period 
of 0.132919(16) d (Fig. 5a) showed a shght phase offset re- 
quiring a period increase of 0.000025 d (— 1.6(t using original 
errors). Plots of the GO-8219 data using the original period 
(Fig. 5c) showed that after ~401 d the phase error had grown 
to nearly 0.6, and by the SM3/ACS-9028 data 1002 days later, 
the phase error had grown to 1 .4. Fortunately, the combination 
of the GO-8267, GO-7503, GO-8219 and SM3/ACS-9028 data 
allows all four data sets to be phased with no cycle count ambi- 
guity. The resulting period and phase when fitting all four sets 
of data are Pb = 0.132944(1) d and T,„^ = 51366.8705(3) MJD 



(see Table 2). By using the archival data the error in the period 
has been reduced by a factor of 16, with no error change in T„ax- 
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Fig. 6. — The difference between the V and / time series from GO-8267 
plotted as a fimction of phase for W29opt. The variable is hottest when it is 
brightest (</ij, = 0.0) and coldest at 4>i, ~0.4. 




70.6 70.8 

Time (HJD 

Fig. 7. — Segment of W34opt time series for the GO-8267 data, showing 
both the V and / time series. A sinusoidal fit to the time series is shown. Note 
the large deviations from this fit near HJD ~51367.8 and HJD = 51370.35- 
51370.5. 

Using these improved results, the periods of 47 Tuc W 
and W29opt now differ by 0.000006 d (0.5 s), a difference of 
0.2 a when using the larger (radio) error. Projecting r„,ai for 
W29opt backwards in time by 1147 orbits, we derive T^ax = 
51214.3837(12) MJD, which differs from the radio-inferred 
value by only (0.0008 ± 0.0012) d, a 1.2 minute difference 
(0.006 in phase). This combination of excellent binary period 
and phase agreement between W29opt and 47 Tuc W ensures 
that W29opt is the optical counterpart to the MSP. 

2.5.3. Color changes 

We plot in Figure 6 the difference between the V and /-band 
relative intensity plots (see Fig. 4) as a function of orbital phase. 
This figure is equivalent to a V -I color vs phase plot except 
that increasing y-values correspond to bluer colors (unlike with 
magnitudes). Note that (/) W29opt clearly has a bluer color at 
<j)b ~0.0 than at other phases (as expected if the secondary is 
heated by irradiation from a luminous disk or MSP; see § 3), 
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Fig. 8. — Phase plots and power spectra for W34opt. As explained in the text the observed period is close to, but statistically different from, the orbital period. 



and (ii) the minimum is at (pb ~0.4, rather than at (j)b = 0.5. 
The main contribution to this possible light curve asymmetry 

(or excess of V-band signal at (f)h ^0.5) is visible in the V-band 
plot of Figure 4, where several consecutive points lie above the 
sinusoidal model near (f)i, = 0.5. However, given the faintness of 
W29opt at these phases, it is difficult to confidently declare that 
this is a real effect. 

2.6. Short Period Variable W34opt 

From a review of the optical variables we note that one other 
star has properties similar to W29opt. A faint (y=22.3) star with 
blue colors (U-V ^0) shows relatively large ampUtude sinu- 
soidal variations in both V and / and is only 4 "9 away from the 
cluster center The fit to the GO-8267 time series of this star 
gives amplitude = 0.301(7), Pt = 0.0676705(78) d, and T^ax 
= 51366.5573(3) MJD. The Hkely Chandra counterpart, W34 
(Grindlay et al. 2001a), lies 0."19 (2 a) away from the optical 
variable. It has a similar X-ray luminosity to that of W29 (see 
Table 2), and is also a moderately hard X-ray source, though 
marginally softer than W29. 

Two portions of the time series for W34opt are shown in Fig- 
ure 7. The large deviations from the sinusoidal model will be 
discussed in § 3. The phase plots and power spectra for W34opt 
are shown in Figure 8 and the color vs phase plot in Figure 
9. Note from the sinusoidal fit (and Fig. 7) that the period 
of W34opt (97.45 min) is close to the HST orbital period (96.4 
min). However, we are convinced that the variations of W34opt 
are real because: (?) the W34opt and HST orbital periods dif- 
fer at almost the 100 cr level, (/;) in cases where artifacts of 
the HST orbit are present, there is usually significant power at 
the harmonics of the //5r period and at ±(1/86400) Hz, but for 
W34opt the time series are very clean at the one frequency, (/;;) 
since the beat period between 97.45 min and 96.4 min (6.2 d) is 
shorter than the total observing window (8.3 d), there is a full 
phase rotation of one cycle in the W34opt signal with respect to 
the orbit of HST, and {iv) the mean difference image (difference 



image at (j)b ~0.0 minus the difference image at ~0.5) has 
a stellar shaped PSF, which is only expected if the variation is 
real, rather than an artifact from a bad pixel interacting with the 
HST orbit-induced motions to generate a false signal. 




Binary phase 

Fig. 9. — The difference between the V and / time series (from GO-8267) 
plotted as a function of phase for W34opt. 



3. DISCUSSION 



3.1. W29opt 

As shown above, there is excellent period and phase agree- 
ment between W29opt and 47 Tuc W and we regard the asso- 
ciation between these two objects as secure. This marks the 
first time a radio pulsar has been localized with sub-arcsecond 
precision by means other than radio pulse timing or synthesis 
observations. We now discuss the nature of W29opt and its vari- 
ability in the context of all that is known about the binary sys- 
tem, including radio (§ 2.1) and X-ray data. 

As noted in § 2.1, the MSP 47 Tuc W was eclipsed for ~30% 
of its orbit on the only day in which it was detected at radio 
wavelengths. The two eclipsing MSPs known in the field (PSRs 
B1957-H20 and J205 1-0827, with companion masses of 0.02- 
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0.03 M0) show, to first order, persistent eclipses at similar or- 
bital phases (other ecUpsing MSPs with similar binary parame- 
ters, in 47 Tuc, are far more poorly studied and relatively little 
is known in detail about their eclipsing properties). By contrast, 
the MSPs in NGC 6397 (NGC 6397 A; D'Amico et al. 2001) 
and Terzan 5 (PSR B1744-24A; Nice & Thorsett 1992) show 
eclipses that are variable in time and orbital phase. These two 
systems, while significantly different from each other, are in ad- 
dition distinct from the other eclipsing MSPs in having larger 
companion masses, ^0.45 M0 and ^0.10 Mq, respectively. 

Naturally we do not know if the eclipses in 47 Tuc W are 
persistent, but the minimum companion mass of 0.13 M0 sug- 
gests it is a fundamentally different star from the ^^0.02 Mq 
objects. The companion mass is possibly similar to that of 47 
Tuc U's companion and other presumed He WD companions 
to MSPs, but the observed eclipse could not be caused by a He 
WD companion, and in any case as seen in Figure 2 the pho- 
tometry for W29opt appears inconsistent with a He WD. The 
companion mass is above the hydrogen burning limit. It there- 
fore appears plausible that W29opt is an MS star, the first such 
known example as a companion to an MSP (if the companion 
in NGC 6397 A is an evolved subgiant rather than an MS star). 

The large amplitude, nearly sinusoidal light curve of W29opt 
is reminiscent of the optical fight curves of PSRs B1957H-20 
(Callanan et al. 1995) and J205 1-0827 (Stappers et al. 2001). 
This similarity suggests that the optical variations of W29opt 
are caused by rotational modulation, as we see alternately the 
heated and (relatively?) unheated sides of the tidally locked 
companion once per orbit (as expected, W29opt appears to be 
hottest when it is brightest; Fig. 6). A similar, large am- 
plitude light curve may also be present in the companion to 
PSR B 1744-24A in Terzan 5, but the formidable reddening and 
crowding in this cluster (Cohn et al. 2002) have thwarted a re- 
cent attempt to detect the companion with HST (Edmonds et al. 
2001b). 

Since the blue color of W29opt probably results from the ef- 
fects of heating, we have estimated the effects and intensity of 
the 'radiation bath' W29opt experiences, based on the size of the 
binary and the expected spin-down luminosity (E) of 47 Tuc W. 
Unfortunately, without a timing solution, no direct estimate of 
E is possible. Using the observed X-ray luminosity of 7.8 x 10^° 
erg s"' and the Lx vs E relationship of Grindlay et al. (2002), 
we estimate that E for 47 Tuc W is ~ 3 x 10^^ erg s'^ (the 1 a 
range is ^ 2-8 x lO"'^ erg s"'). A much smaller value of E 
(7.8 X lO-'^ erg s~^) is derived using the linear relation of Becker 
& Trumper (1999), L^r^lQ-^E, derived largely for field MSPs. 

The stellar evolution models of Bergbusch & Vandenberg 
(1992) give an approximate guide to the relationship between 
luminosity and mass for a MS star We adopt the lowest mass 
model given by Bergbusch & Vandenberg (1992), a 0.15 Mq 
model with Tgff = 3300 K and My = 13.2, and test whether sig- 
nificant heating of such a star is expected. We use the mea- 
sured orbital period, an assumed mass for the NS of 1.4 Mq, 
and the above secondary mass to estimate the binary separation 
from Kepler's Third Law. Then, using the estimated secondary 
radius from the stellar models, and assuming that the MSP ra- 
diates its wind isotropically, the 'high' E estimate given above 
implies that the MSP energy intercepted by W29opt should be 
~ 1 X lO^^erg s~', a factor of ~ 120x the luminosity of the as- 
sumed companion (by comparison, the corresponding numbers 
for the companion to 47 Tuc U are 7 x lO'^'erg s"^ or 0.013 x 
the luminosity of the companion; EGHOl). Assuming that aU 



of the spin-down energy is reradiated as a blackbody, we de- 
rive a temperature of 15,500 K. This would imply that heating 
effects are important and the MS star should be brighter and 
hotter than in its unperturbed state, possibly explaining the blue 
colors and relative brightness of W29opt. 

The above temperature is an upper limit (unless there is 
beamed emission), since a significant fraction of E could be 
powering other processes such as mass loss. For example, for 
PSR 1957H-20, CaUanan et al. (1995) find that the fraction of 
the MSP flux intercepted by the secondary that is converted into 
optical emission (defined as rj) is either in the range 0.07-0.2 
for a Roche lobe filling model, or > 3 (significant beaming) 
for a secondary considerably underfilling its Roche lobe. Stap- 
pers et al. (2001) find t] ~0.3 for J205 1-0827 and EGHOl find 
r; ~ 1 for 47 Tuc U. Guided by these numbers we assume a con- 
servatively small value for r] of ~0.1 to give a reasonable lower 
limit on the heating flux, and derive a blackbody temperature of 
^8,700 K. Therefore heating effects should still be important. 

For this model, the ratio between the heating flux and the 
undisturbed luminosity of the assumed stellar companion (120) 
is only a factor of 2.4 larger than the ratio (caused by heating) 
between the V-band flux of W29opt and the V-band flux of the 
stellar model (implying tj ^0.4). However, if this model is ap- 
proximately correct, the relatively small difference in intensity 
between the heated and unheated sides of W29opt (a factor of 6 
compared to the above factor of 120) has to be explained. One 
possibility is that heating of the companion has occurred fairly 
evenly over its surface. Although the secondary is probably 
tidally locked, significant heating of the darker side of the star 
might be expected because of energy transport within the star 
from the combined effects of convection (this low mass MS star 
should be fully convective), and the effects of relativistic e"'"e~ 
pairs and 7-ray photons (D'Antona 1996). It is also possible 
that a crescent of the heated hemisphere is visible at all orbital 
phases if the orbital inclination of the binary is significantly 
less than 60° (Fruchter et al. 1995). This requires a higher 
mass companion for consistency with the radio results and the 
brighter star also gives a correspondingly smaller amplitude. In 
the case of NGC 6397 A, where the inclination is known to be 
high, the companion also has a much smaller variability am- 
plitude than expected based on heating estimates (Orosz & van 
Kerkwijk 2002). This may be a feature common to all MSP 
systems containing non-degenerate companions. 

Approximate upper fimits to the mass and radius of the 
47 Tuc W companion are set by the requirement that mass 
transfer be avoided. Using the Roche-lobe formula from 
Paczyhski (1971) (r/a = 0.462[(Mw/(Mns+Mw)]'/\ where r 
is the Roche-lobe radius, a is the binary separation, Mw is 
the mass of the secondary and Mns is the mass of the NS), 
we derive that, in its unperturbed state, a 0.29 Mq model with 
Teff= 3680 K and My = 10.8 would underfill its Roche lobe by 
~1% (the corresponding upper Umit on the orbital incUnation 
is 27°). Our 0.15 Mq model would (in its unperturbed state) 
significantly underfill its Roche lobe (the ratio of the stellar ra- 
dius to r is 0.6), explaining why the NS is not accreting. How- 
ever, under the influence of the extreme radiation bath from the 
MSP, W29opt may expand to fill or almost fill its Roche lobe 
(D'Antona & Ergma 1993), with accretion being inhibited by 
radiation pressure from the MSP (as may be occuring in NGC 
6397 A; Burderi et al. 2002). 

Using the 'low' E estimate given above, the MSP energy in- 
tercepted by W29opt is only 4 x 10^'erg s"' (with a reradiated 
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blackbody temperature of 6,300 K), a factor of ~3.3 times the 
luminosity of the assumed companion. This luminosity ratio is 
smaller than the observed variation and is a factor of 15 smaller 
than the ratio between the V-band flux of W29opt and the V- 
band flux of the stellar model. So, unless significant beam- 
ing is present, only a fainter MS star than assumed above can 
give a larger ratio between the heating flux and the stellar lumi- 
nosity. However, the disparity between the V-band fluxes then 
increases by the same amount. A smaller ratio between the V- 
band fluxes is given by a brighter, more massive secondary (re- 
quiring an incUnation < 60°), but then the heating flux would 
be too small to cause the observed variations. These inconsis- 
tencies suggest that the higher E value given earlier is closer to 
the true value, providing indirect support for the Lx vs E rela- 
tionship of Grindlay et al. (2002). 

The similarities between the X-ray properties of W29 and 
those of the X-ray counterparts of NGC 6397 A (see § 2.3) and 
47 Tuc J offer valuable insight. These similarities at the very 
least allow for a similar physical origin of the X-ray emission. 
The extended nature of the hard X-ray source in NGC 6397 A 
suggests it is due to shocked gas Ufted from the binary compan- 
ion (Grindlay et al. 2002), and similar behavior may be occur- 
ring in W29. This match in X-ray properties between W29 and 
two other echpsing MSPs suggests evidence for a trend, since 
apart from 47 Tuc J, aU 9 of the X-ray-detected MSPs in 47 Tuc 
that are bright enough to have useful hardness ratios determined 
have soft spectra (caused by likely thermal emission from their 
polar caps; Grindlay et al. 2002). 

As we have already mentioned, 47 Tuc W is an extremely 
unusual (and so far possibly unique) MSP in having a Ukely 
MS companion. This likely represents evidence for stellar in- 
teractions, not surprising for a system so close to the cluster 
center. One possible formation scenario for this system is that 
an NS was spun up by accretion to form an MSP, which even- 
tually evolved into an MSP/very low mass companion binary. 
Later, a collision between the MSP binary and the currently ob- 
served MS star caused the very low mass star (the lightest object 
amongst the 3 stars) to be ejected from the system. Alterna- 
tively, a single MSP could have exchanged into a double MS 
star binary, although in this case tight constraints on the mass 
of the ejected star may be necessary unless W29opt was once 
considerably more massive than it is today. Significant mass 
loss could be taking place if W29opt is close to fiUing its Roche 
lobe and material is being expelled from the system under the 
propeller mechanism (Burderi et al. 2002). A third possibility 
is that the observed binary was formed as a result of direct 2- 
body tidal capture (Mardling 1995) between a single MSP and 
an MS star, although this process may be much less important 
than binary interactions (Rasio et al. 2000). 

In testing these various exchange scenarios, we note that 47 
Tuc W is located very close to the cluster center (only 3 "8 or 
0.16rc away), consistent with the centrally concentrated spa- 
tial distribution of the other MSPs. If the MSP was initially 
ejected from the cluster core in an exchange encounter it may 
have had time to sink back into the cluster core, since the MS 
star could have been captured several Gyr ago, and the cen- 
tral relaxation time for 47 Tuc is 0. 1 Gyr (increasing to 3 Gyr 
at flie half-mass radius of 2.'8; Harris 1996). The NGC 6397 
MSP, by contrast, is found ^30" from the cluster center (Fer- 
raro et al. 2001) and is therefore well outside the cluster core 
(fc = 3."0; Harris 1996), possibly showing evidence for ejection 
from the core (D'Amico et al. 2001; Grindlay et al. 2002). 



This ejection may have occurred quite recently, since the cen- 
tral relaxation time for the core-collapsed NGC 6397 is only 
0.08 Myr, increasing to 0.3 Gyr at the half-mass radius of 2.'3 
(Harris 1996). If these numbers are accurate, a system as heavy 
as NGC 6397 A should sink rapidly back into the core (unless 
it is in a highly elliptical orbit). One possible explanation is that 
NGC 6397 A was formed very recently by an NS capturing a 
subgiant, and was ejected from the core in this interaction. This 
scenario is broadly consistent with the relatively small (charac- 
teristic) age for the system of 0.35 Gyr (D'Amico et al. 2001). 
It may also have been ejected from the cluster core in a sub- 
sequent interaction. A similar ejection scenario may apply to 
PSR B1744-24A, which is located - 40" from the center of 
Terzan 5 (Nice & Thorsett 1992), equivalent to ^ 5rc using the 
recent determination by Cohn et al. (2002). 

It may be possible that W29opt will eventually fill its Roche 
lobe and cause mass accretion onto the NS to recommence, as 
may occur for NGC 6397 A (Ferraro et al. 2001 and Burderi 
et al. 2002). In this scenario, as W29opt keeps losing mass, it 
eventuaUy may end up like the very low mass (0.02-0.03 Mq) 
systems such as 47 Tuc J (though these may be ablated He WDs 
rather than ablated MS stars; Rasio et al. 2000). Then, after 
continued ablation, it could possibly become an isolated MSP. 



Unlike with W29opt, there is no period or phase match for 
W34opt with a known MSP (its period is similar to, but statis- 
tically different from, the 95.39 ±0.1 min period of 47 Tuc R, 
an echpsing MSP with a very low mass companion; CLFOO), 
but a significant fraction of the 47 Tuc MSPs have not yet been 
detected in the radio (CLFOO; Grindlay et al. 2002). Given 
the short period of this system, the heating effects (if it is an 
MSP) may be even more dramatic than for W29opt, since the 
ratio between the heating flux and the stellar luminosity would 
be higher, assuming a similar luminosity companion star. We 
note that the color difference between the brighter and darker 
sides of the companion (^30%) compared to the amplitude of 
the variation is larger than the similar ratio for W29opt. The 
smaller amplitude of W34opt compared to W29opt may simply 
be an inclination effect or it may reflect a significant difference 
in the level of heating of the dark side of the companion, or 
show evidence for a different type of star, such as a very low 
mass degenerate star instead of an MS star. We also note that 
in this system the low mass model assumed above for W29opt 
would overfill its Roche lobe by a few percent. 

There are broad similarities in the appearance of the color vs 
phase plots for W29opt (Fig. 6) and W34opt (Fig. 9), particu- 
larly in the positions of the peak near (pi, = 0.0, the minimum 
at (ph '^OA, and the possible 'hump' of relatively blue light at 
(pb = 0.4-0.7. The relatively blue color at optical maximum 
is secure, but the effects near (ph = 0.5 can only be considered 
marginal, given the faintness of the stars at these times and the 
possibiUty that artifacts of the HST orbit may have leaked into 
the time series (half the period of W29opt is 95.72 min). 

These comparisons do not constitute proof that W34opt is an 
MSP companion, and therefore we briefly consider alternative 
explanations for its behavior. The near-sinusoidal Ught curve 
of W34opt, with its short period and large amplitude, is sim- 
ilar to that of a W UMa variable (a contact, or near contact 
binary consisting of two MS stars). However, this interpreta- 
tion is clearly inconsistent with the blue colors of this object. 
A second possible explanation is that W34opt is a CV, since 
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most of the CVs in 47 Tuc (and also most of the CVs in NGC 
6397 and the field) are optically variable blue stars with moder- 
ately hard X-ray counterparts (Grindlayetal. 2001a). However, 
the likely CVs in 47 Tuc show flickering (random fluctuations 
on timescales of minutes with amplitudes of '^0.05-0. 1 mag) 
which hides the presence of periodic variations in most cases. 
Of the few CVs in 47 Tuc which do show periodic variations, 
the light curves are very different from those of W34opt (and 
W29opt), and show either low amplitude (< 0. 1 mag) periodic 
variations from ellipsoidal modulation, or a combination of el- 
lipsoidal modulation and eclipses (Edmonds et al. 2002a, in 
preparation), with distinctly non-sinusoidal light curves (there 
is significant signal in the power at both the orbital period and 
half the orbital period). Given these differences, then unless 
W34opt is a different type of CV from the ^^20 good candidates 
found in 47 Tuc, we consider it more likely to be an MSP. 

As shown in Figure 7, there are several instances of large 
deviations from the sinusoidal model for W34opt. We exam- 
ined the standard deviation of the residuals (data-sinusoid) for 
each cycle of the W34opt light curve, after iteratively remov- 
ing 3 a deviations from the mean, and compared these results 
with those found for W29opt. The mean standard deviation for 
W34opt (0.15) was similar to the value found for W29opt (0.21), 
but there were some differences in the distributions. No 3 a de- 
viations were found for W29opt, but two were found for W34opt 
(at HJD ~51363.6 and HJD ~51367.8), and two consecutive 
cycles of ~ 2.9a deviations (in HJD = 5 1 370.35-5 1370.5) were 
also found. These variations could be evidence for variable 
mass loss from the companion or they could be episodes of un- 
usually large flickering, if the CV model is correct. 



4. SUMMARY AND FUTURE OBSERVATIONS 

Briefly summarizing the key results in this paper, we have 
identified a binary companion, W29opt, to a second MSP in the 
cluster 47 Tuc (adding to the detection of a He WD companion 
to MSP 47 Tuc U; EGHOl). This faint, blue star shows large 
amplitude sinusoidal variations in both V and / (caused by heat- 
ing of a tidally locked companion), and has a period and phase 
which match remarkably well with those of the MSP 47 Tuc W. 
We argue, based on the radio and optical results, that W29opt is 
a main sequence star, a rare circumstance for an MSP compan- 
ion, and we present evidence for another optical variable (with 
similar properties to W29opt) that may also be an MSP compan- 
ion. 

Deeper Chandra observations of 47 Tuc scheduled for 2002 
September (PI J. Grindlay) will allow variability and spectral 
studies of W29, with dependences on phase, to be made. Simul- 
taneous observations scheduled with ACS (as part of the same 
program) wiU provide the first EST images of 47 Tuc in Ha, 
and should detect evidence for a strong emission line in W34opt 
if it is a CV rather than an MSP companion. Combination of 
this data with archival data and with the GO-8267 time series 
will refine the period and phase of W34opt. Finally, by combin- 
ing these X-ray data with existing and future //-ST observations 
of 47 Tuc, more MSP companions may be identified. 

We thank Andy Fruchter and Kailash Sahu for discussions, 
and the referee Scott Ransom for helpful comments. This work 
was supported in part by STScl grant GO-8267.01-97A (PDE 
and RLG) and NASA grants NAG 5-9095 and SAO grant GOl- 
2063X (EC). 
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An MSP Optical Counterpart in 47 Tuc 



Table 1 
Summary of optical data 



Dataset 


Observation Date 


Camera 


Filters 


Exposures 


GO-8267 


1999 July 3-11 


WFPC2 


F555W CV) 


636 x160 s 




WFPC2 


F814W ('/') 


653 x160 s 


GO-7503 


1999 October 28 


WFPC2 


F555W 


9x40 s; 3x120 8" 


GO-8219 


1999 September 10 


STIS 


CLEAR 


8x30 s 




2000 August 7-16 


STIS 


CLEAR 


38x30 s 


SM3/ACS-9028 


2002 April 5 


ACS/HRC 


F475W 


20x120 s 



'The 40 s exposures were composed of co-added 20 s exposures 



Table 2 

Optical, radio and X-ray data for the companion to MSP 47 Tuc W (W29opt) and a candidate MSP counterpart (W34, 



Source 


(J2000) 


Dec'' 
(J2000) 


U 


V 


/ 


(days) 


(MJD) 


(ergs-l) 


X-ray spectrum'^ 
(photon index) 


W29op, 


00 24 06.07(1) 


-72 04 49.02(6) 


23.8(5) 


22.3(2) 


21.6(3) 


0.132944<1)*^ 


51366.8705(3) 


7.8X10™ 


1.8±0.6 














51214.3837(12)^ 






47 Tuc W 












0.13295(4)f 


51214.38454(3)8 


4X10™ 


2 4+0.5 


W34op, 


00 24 05.21(2) 


-72 04 4 659(6) 


22.3(5) 


22.3(3) 


21.6(4) 


0.0676705(78) 


513665573(3) 



^Coordinates are in the JPL DE200 planetary ephemeris frame used for MSP timing (Freire etal. 200ia) 
''Time of optical maximum (see text) 

^Lx and die spectral fit assume cluster iVfl = 3.0 X lO^^cm"^ 

"^Binary period derived &om optical data 

^Oplical Tmax - 1 147 orbits 

^Binary period derived from radio data 

S Tfffax derived from radio data 



